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ABSTRACT 

We investigate the influence of blending on the Cepheid distance scale. Blending 
leads to systematically low distances to galaxies observed with HST, and therefore 
to systematically high estimates of the Hubble constant Hq. We select a sample of 
43 long-period, large-amplitude Cepheids in the Large Magellanic Cloud, from the 
catalog recently released by the Optical Gravitational Lensing Experiment. We then 
model the effects of blending, as observed by HST at large distances, by adding the 
V, /-bands contributions from nearby bright stars. We find that the derived distance 
modulus would be too short, compared to the true value, by ~ 0.07 mag at distance 
of 12.5 Mpc and by ~ 0.35 mag at distance of 25 Mpc. This has direct and important 
implications for the Cepheid distances to galaxies observed by the HST Key Project 
on the Extragalactic Distance Scale and other teams: half of the KP sample is likely 
to exhibit a blending bias greater than 0.1 mag. 



Introduction 



As the number of extragalactic Cepheids discovered with HST continues to increase and 
the value of Hq is sought from distances based on these variables (Gibson et al. 1999, Saha et 
al. 1999), it becomes even more important to understand various possible systematic errors which 
could affect the extragalactic distance scale. Currently, the most important systematic is a bias in 
the distance to the Large Magellanic Cloud (LMC), which provides the zero-point calibration for 
the Cepheid distance scale. The LMC distance is very likely ~15% shorter than usually assumed 
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(e.g. Udalski et al. 1999a; Stanek et al. 1999), but it still might be considered uncertain at the 
~ 10% level (e.g. Jha et al. 1999). Another possible systematic, the metallicity dependence of 
the Cepheid Period-Luminosity (PL) relation, is also very much an open issue, with the empirical 
determinations ranging from to —0.4 mag dex~^ (Freedman & Madore 1990; Sasselov et al. 1997; 
Kochanek 1997; Kennicutt et al. 1998, Udalski et al. 1999a). 

In this paper we investigate a much neglected systematic, that of the influence of blended 
stellar images on the derived Cepheid distances. Although Cepheids are very bright. My ~ —4 at 
a period of 10 days, their images when viewed in distant galaxies are likely to be blended with 
other nearby, relatively bright stars. Recently Mochejska et al. (1999) showed that a significant 
fraction of Cepheids discovered in M31 by the DIRECT project (e.g. Kaluzny et al. 1998; Stanek 
et al. 1998) were resolved into several blended stars when viewed on the HST images. The average 
FWHM on the DIRECT project ground-based images of M31 is about 1.5", or ~ 5 pc, which 
corresponds to the HST-WFPC2 resolution of 0.1" for a galaxy at a distance of 10 Mpc. Any 
luminous star (or several of them) in a volume of that cross section through a galaxy could be 
indistinguishable from the Cepheid and would contribute to its measured flux. 

In this paper we investigate the effects of stellar blending on the Cepheid distance scale by 
studying Cepheids and their close neighbors observed in the LMC by the Optical Gravitational 
Lensing Experiment (OGLE: Udalski, Kubiak & Szymahski 1997). The catalog of ~ 1300 LMC 
Cepheids has been recently released by Udalski et al. (1999b). As the LMC is ~ 100 — 500 times 
closer to us than galaxies observed by HST, ground-based resolution of ~ 1.0" allows us to probe 
linear scales as small as ~ 0.25 pc in that galaxy. 

We describe the OGLE data used in this paper in Section 2. In Section 3 we apply these data 
to simulate the blending of Cepheids at various distances. In Section 4 we discuss the implications 
of our results for the Cepheid distance scale. In Section 5 we propose further possible studies to 
learn more about the effects of blending. 



2. The Data 

The data used in this paper came from two catalogs produced by the OGLE project. The 
first one, with 1333 Cepheids detected in the 4.5 square degree area of central parts of the LMC, 
has just been released (Udalski et al. 1999b) and it is available from the OGLE Internet archive 
at http://www.astrouw.edu.pl/~ftp/ogle/. It contains about 3.4 x 10^ BVI photometric 
measurements for the variables, along with their derived periods, BVI photometry, astrometry 
and classification. 

The second catalog, that of BVI photometry for many millions of LMC stars observed by the 
OGLE project, will be released soon (Udalski et al. 1999c, in preparation). Its construction will 
be analogous to the SMC BVI maps (Udalski et al. 1998). A typical accuracy of the photometric 
zero points of the LMC photometric data is about 0.01 — 0.02 mag for all BVI-hands and the 
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catalog reaches V ^ 21 mag. 

As our goal was to estimate the influence of blending on the Cepheid distance scale, we 
selected for further investigations only 54 longest period {P > 10 days), fundamental-mode 
Cepheids. We further required their /-band total amplitude of variations to exceed 0.4 mag, 
which corresponds to amplitude of ~ 0.7 mag in the l^-band. This is to reflect the fact that in 
distant galaxies a typical photometric error of ~ 0.1 mag prevents discovery of low amplitude 
Cepheids. We used the /-band amplitude criterion because OGLE observes mostly in the /-band. 
The amplitude cutoff reduced the OGLE LMC sample to 47 Cepheids. Finally, we have excluded 
four highly reddened Cepheids, which left us with 43 variables. It should be noted that because 
of the CCD saturation limits in the OGLE data the longest period Cepheid in our sample has 
P = 31 days. This sample will be used to investigate the effects of blending on the Cepheid 
distance scale. 



3. Cepheids and Their Neighbors 

The LMC is located at about 50 kpc, or fJ-o^LMC = 18.5, from us (for simplicity, in the rest 
of this paper we use the distance scale as adopted by the HST Key Project on the Extragalactic 
Distance Scale, hereafter: KP). We define two distances for which we simulate the blending using 
the LMC data: 12.5 Mpc {fiQ ~ 30.5), somewhat shorter than the median distance of 14.4 Mpc 
for the KP sample (see Table 1 in Ferrarese et al. 1999), and 25 Mpc {fiQ ^ 32.0), roughly 
corresponding to the most distant galaxies in which Cepheids were detected with HST (Gibson et 
al. 1999). The FWHM of HST-WFFC2 is ~ 0.1", which at 12.5 Mpc and 25 Mpc corresponds 
physically to area of radius 12.5" and 25" at the LMC distance. 

In Figure || we show two OGLE LMC Cepheids and their neighbors in the /-band (left 
panels) and the l^-band (right panels). The images are 1' in size. Also shown are two circles 
corresponding to the FWHM of HST-WFFC2 camera at 12.5 Mpc and 25 Mpc, i.e. 12.5" and 25" 
in radius. The two Cepheids were chosen to represent two different situations: the one shown in 
the top panels, LMC_SC15 118594, has only one bright and nearby, red neighbor at 5.4" from the 
Cepheid, and several other, fainter neighbors further away. The second one shown in the bottom 
panels, LMC_SC11 250872, is located in a very dense stellar region and it is probably a member 
of a small star cluster. Part of the cluster light would be included in HST-WFPC2 measurements 
of the Cepheid if the LMC were at 12.5 Mpc while at 25 Mpc almost entire cluster light would 
be included. As discussed later in the paper, large amount of blended light would most probably 
cause LMC_SC11 250872 to elude detection when observed at large distances. 

The OGLE photometric catalog of stars in the LMC extends about six magnitudes, or -^0.4% 
in flux (see Figure 4 of Udalski et al. 1999b) below the bright sample of Cepheids selected in the 
previous Section. We want to define a criterion to separate stars which will contribute to the flux 
of a Cepheid when blended together, from those which would contribute only to the background 
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Fig. 1. — Two LMC Cepheids and their neighbors shown in the /-band (left panels) and the V- 
band (right panels). The images are l' in size. Also shown are two circles corresponding to the 0.1" 
FWHM of HST-WFPC2 camera at 12.5 Mpc and 25 Mpc, translating into radius of 12.5" and 25" 
at the LMC distance. One pixel of the OGLE camera is 0.417". 

light in the host galaxy. We use a lower limit of 5% of the mean flux of the Cepheid for a star to 
be included as a blend, but will discuss different values later in the paper. 

We use this 5% cutoff in evaluating the sum Sp (Mochejska et al. 1999) of all flux contributions 
in filter F normalized to the flux of the Cepheid: 

Np , 
i=l Jc 

where is the flux of the i-th companion, fc the mean intensity flux of the OGLE LMC Cepheid 
and Np the total number of companions within either 12.5" and 25" in radius. In Figure |2| we 
show the cumulative probability distribution of flux contribution from companions Si (left panels) 
and Sy (right panels) within a radius of 12.5" and 25" of LMC Cepheids. For the smaller radius 
of 12.5" 20-25% of our sample is not blended and the contribution of blue blends is somewhat 
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Fig. 2. — The cumulative probability distribution of flux contribution from companions Si (left 
panels) and Sy (right panels) within a radius of 12.5 and 25" of LMC Cepheids, corresponding to 
0.1" FWHM of HST-WFPC2 camera at 12.5 (upper panels) and 25 Mpc (lower panels). 5/,5y 
are defined (Equation 1) so that when Sj, Sy = 1, the total flux from the blended neighbors equals 
to the mean flux from the Cepheid. 

stronger than that of red blends. For the larger radius of 25" all 43 Cepheids are blended to some 
extent and the contribution of red blends is now more significant than that of blue blends. In 
the next Section we attempt to use our data to quantify the effects of blending on the Cepheid 
distance scale. 



4. Blending and the Cepheid Distance Scale 

Wc adopt the KP procedure for deriving distances, as given by prescription in Madore & 
Freedman (1991). LMC is assumed to be at distance modulus of l-to,LMC = 18.50, with LMC 
Cepheids reddened by E{B — V) = 0.10 mag. When we apply this prescription to the V,I data 
of our 43 OGLE Cepheids, we obtain values of ij.o,lmc = 18.56 and E{B — V)= 0.12 mag, i.e. 
somewhat discrepant, but basically indicating fairly good agreement in photometry. The rms 
scatter around the fitted P-L relations (with their slopes fixed to that of the KP prescription) is 
0.17 mag in the F-band and 0.13 mag in the /-band. 

As the next step we add contributions from the nearby stars to each Cepheid (applying the 
5% cutoff discussed in the previous section) in the V and I band separately and we repeat the 
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Fig. 3. — Blending effects simulated at 12.5 (left panels) and 25 Mpc (right panels). The continuous 
lines show the best-fit PL relations when the blends are added, and the dashed lines those without 
the blending. For discussion see text. 



distance derivation procedure. To take into account the fact that a heavily blended Cepheids 
would elude detection, we require that their blended /-band total amplitude of variations exceeds 
0.4 mag. The results for simulated distances of 12.5 and 25 Mpc are shown in Figure ^. This 
Figure shows a number of interesting features and deserves detailed discussion. 

For the simulated distance of 12.5 Mpc, the sample of Cepheids is reduced to 35 and the 
derived distance modulus is smaller than the "true" (unblended) value of Hq^lmc = 18.56 by 
0.07 mag. The reddening estimate is E{B — V) = 0.10 mag, smaller than for the unblended sample 
because of the slightly larger contribution of blue blends, as discussed in the previous Section. 
The rms scatter around the fitted P-L relations increases to 0.22 mag in the l^-band and 0.17 mag 
in the I-band. This is because while there are now Cepheids in the sample with fairly substantial 
blending, it still contains Cepheids with no blending (see Figure |^ . 

The situation becomes quite dramatic at the simulated distance of 25 Mpc (right panel 
of Figure ^). There are only 13 Cepheids left with /-band amplitude larger than 0.4 mag, 
with the shorter period (and therefore fainter) Cepheids preferentially removed. The derived 
distance modulus is smaller than the "true" value by 0.36 mag, with the reddening estimate 
E{B — y) = 0.11 mag. What is very interesting is that the rms scatter around the fitted P-L 
relations now decreases to 0.17 mag in the y-band and 0.13 mag in the /-band. This is because 
the sample of Cepheids, while much smaller, is now more homogeneous when it comes to blending. 



0.4 



0.3 



CO 



^ 0.2 



3, 



0.1 



■ blends > 10% 
O blends > 5% 
A blends > 2.5% 



~i — I — I — r 



m 

O 



O 



A 
O 



o 



A 



25 



5 10 15 20 

d [Mpc] 

Fig. 4. — Blending difference between the true and measured distances modulus shown as a function 
of simulated distance. Three different cutoffs for blending are investigated (see text for discussion): 
10, 5 and 2.5% of the mean intensity flux of the Cepheid. Also shown is the histogram of distances 
for 24 galaxies for which Cepheid i/S'T distances were measured or re-measured by the KP (Ferrarese 
et al. 1999). 



All Cepheids are to some extent blended (see Figure with strongly blended cases removed by 
the high amplitude requirement. 

Since there is such a dramatic difference in blending between the two distances simulated 
so far, we decided to investigate the blending for a larger number of distances. The results are 
shown in Figure ^, where blending difference between the true and measured distance modulus is 
shown as a function of simulated distance. Also shown is the histogram of distances for 24 galaxies 
for which Cepheid distances with HST were measured or re-measured by the KP (Ferrarese et 
al. 1999). 

The 5% cutoff which we employed to define blended stars is somewhat arbitrary and in reality 
is most likely a function of data reduction procedure, signal-to-noise in the images etc. We decided 
to investigate the dependence of blending effects on the cutoff value by using two additional 
cutoffs: 2.5% and 10%. The results are shown with different symbols in Figure H Clearly, the 
exact value of the blending difference between the true and measured distances modulus depends 
on the applied cutoff, but the overall trend remains the same. 

Looking at the histogram of distances in Figure ^ we can see that half of the KP sample is 
likely to exhibit a blending bias greater than 0.1 mag, and in some cases it can be as large as 
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0.3 mag. Clearly blending influence on the Cepheid distance scale can be potentially very large 
and cannot be neglected. 

5. Further Studies of Blending and Conclusions 

The study of Mochejska et al. (1999) showed that a significant fraction of Cepheids discovered 
in M31 by the DIRECT project were resolved into several blended stars when viewed on the 
HST images. As we have shown in this paper, modelling of the blending effects using Cepheids 
in the LMC possibly indicates a large, 0.1 — 0.3 mag bias when deriving Cepheid distances to 
galaxies observed with HST. In addition, blending is a factor which always contributes in only one 
direction, and therefore it will not average out when a large sample of galaxies is considered. The 
sign of the blending effect on the Hq is opposite to that caused by the lower LMC distance (e.g. 
Udalski et al. 1999a; Stanek et al. 1999). 

We would like to point out that the blending of Cepheids is likely not only to affect the 
studies of these stars in different galaxies, but might also affect differential studies, such as that 
of Kennicutt et al. (1998) in the spiral galaxy MlOl. Their observed effect that metal-rich (and 
therefore closer to the center of the galaxy) Cepheids appear brighter and closer than metal-poor 
(and therefore further away from the center) stars could be partially caused by the increased 
blending closer to the center of the galaxy, although at this point we make no attempt to estimate 
how much of this effect would be indeed due to blending. 

The bar of the LMC, where most of data discussed in this paper were collected by OGLE, 
seems on average to have higher surface brightness than a typical KP galaxy (Macri et al. 1999). 
It would be desirable to further establish the importance of blending for the Cepheid distance 
scale using a variety of methods and data in a number of different galaxies. Mochejska et al. (2000, 
in preparation) are now studying the HST archival images of a large sample of ~ 100 Cepheids 
detected in M33 (Macri et al. 2000, in preparation) by the DIRECT project. An approach 
analogous to that used in this paper will be employed by Stanek & Udalski (1999, in preparation) 
for a sample of OGLE Cepheids in the Small Magellanic Cloud, which will have the advantage of 
including Cepheids with periods of up to P ~ 50 days in a lower surface brightness system (de 
Vaucouleurs 1957). 

Another approach, closer reproducing the procedure employed by the KP when using 
Cepheids to determine distances, would be to use HST images of relatively nearby galaxies, such as 
NGC3031 or NGC5253 (Ferrarese et al. 1999), by degrading them in resolution and signal-to-noise 
as to represent much more distant galaxies. Unfortunately, much of the data for the several closest 
galaxies have been taken before the refurbishment of HST and are therefore of inferior spatial 
resolution compared to later WFPC2 data. 

All these studies can provide only an approximate answer to the blending problem, as 
each individual galaxy can in principle be different in its blending properties. One would like 
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to find a way to constrain or eliminate blending in each individual case. As pointed out by 
Mochejska et al. (1999), one obvious solution to the problem of blending would be to obtain 
data with better spatial resolution, such as planned for the Next Generation Space Telescope 
{NGST). While there will be desire to use NGST to observe much more distant galaxies than 
with HST, it would be of great value to study some of the not-so-distant ones as well. Another 
possible approach would be to try to circumvent the blending problem by developing and 
applying a Period-Amplitude-Luminosity (PAL) relation for Cepheids (Paczyhski 1999, private 
communication), together with image subtraction techniques such as that developed by Alard &: 
Lupton (1998). 
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